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AUTHOR'S  SUMMARY 


The  author  has  chosen  two  aspects,  namely  assessment  by  calculation  and  the 
dynamics  of  composite  materials. 

In  the  first  part,  the  theoretical  studies  are  reviewed,  in  particular 
various  methods  for  evaluating  elastic  constants  (boundary  methods  using 
variational  calculation,  methods  using  Airy  stress  function,  elastodynamic 
methods)  and  recent  methods  of  calculation  for  multilayers  in  the  plastic  state. 
Finally,  a  critical  study  of  the  various  theories  available  for  the  propagation 
of  waves  and  vibrations  in  composite  materials  is  made  (equivalent  homogeneous 
media  theory,  micros tructure  theory,  mixture  theory,  theory  based  on  asymptotic 
development,  etc). 

In  the  second  part,  some  experimental  researches  are  reviewed  for  the 
calculation  of  elastic  constants  in  micropolar  elasticity  under  a  static  and 
dynamic  regime.  The  mechanical  characteristics  of  composite  materials  are 
examined  along  with  the  corresponding  experimental  techniques  (vibration  of  rods, 
ultrasonics,  etc). 
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1  INTRODUCTION 

During  the  past  five  years,  special  composite  materials  -  long  fibres,  in 
particular  -  have  been  developed  for  the  specific  needs  of  the  aeronautical  and 
space  industry.  In  other  sectors  of  industry  more  and  more  use  is  being  made  of 
these  materials  as  complements  to  known  metallic  materials  and  plastics.  These 
products  are  used  each  time  it  is  desired  to  blend  high  mechanical  performance 

with  lightness  of  structure  or  of  mechanical  parts. 

We  are  in  a  phase  where  ability  is  developing  very  fast,  and  often,  in  the 
recent  past,  it  has  outstripped  knowledge.  Around  1960,  this  situation  was  an 
important  factor  for  inciting  researchers.  It  became  transformed,  very  rapidly, 
by  a  veritable  explosion  of  studies  and  researches,  some  theoretical,  others 
experimental,  in  all  fields.  It  is  sufficient  to  peruse  the  technical  reviews  in 
applied  mechanics  to  become  aware  of  this  activity. 

This  interest,  which  borders  on  infatuation,  is  surely  explained  by  the 
desires  of  engineers  and  technicians  to  fill  in  their  basic  knowledge,  which  is 
often  confined  to  conventional  materials.  In  our  opinion  it  is  also  explained  by 
the  fact  that  composite  materials  now  stand  at  a  privileged  crossroads,  to  which 
a  large  number  of  scientific  disciplines  are  converging.  This  situation,  unique 

in  the  history  of  materials,  is  a  very  happy  one  and  quite  rightly  favours  the 

development  of  inter-disciplinary  studies,  as  Fig  1  shows. 

In  fact,  quantum  mechanics,  which  studies  periodic  lattices,  has  studies 

1  2 

which  can  be  transposed  into  the  field  of  composite  materials  ’  .  Thus,  quite 

3 

recently  Kohn  et  at  have  presented  the  results  of  original  studies  in  this  field, 

with  a  variational  formula  only  slightly  known  in  mechanics.  Geology,  which  is 

used  to  treating  geological  layers  in  seismology,  has  studies  which  are  capable 

of  being  extrapolated  to  composite  materials.  In  this  field  are  the  works  of 
4  5  6 

Jardetzky  ,  Brekhovskikh  and  Postma  on  the  propagation  of  waves  in  stratified 
media . 

The  physics  of  metals  and  viscoelastic  materials,  together  with  rheology  in 
general,  are,  naturally,  of  great  interest  in  this  field,  with  the  focus  on  new 
composite  materials,  such  as  orientated  eutecties  and  high  diffusion  fibrous 
materials.  They  can  give  to  the  composite  materials  specialists  a  certain  number 
of  experimental  results  concerning  polyphase  materials  (polycrystalline  metals, 
sequenced  polymers,  amongst  others),  as  Mason  and  Skimin^  show.  In  return, 
recent  studies  on  composite  materials  (periodic  or  random)  could  shed  a  new  light 
on  the  interpretation  of  these  results  (interpretation  ol  multiple  peaks  of 
internal  friction,  wave  attenuation,  etc) . 


PRECEDING?  PAGE  BUkNK-NOT  FlU'tED 


6 


Acoustics  and  optics  have  both  developed,  in  parallel,  numerous  studies  on 
composite  materials.  They  supply  reseaichers  with  numerous  methods  of  investi¬ 
gation,  ranging  from  the  simplest  to  the  most  sophisticated  (for  example,  holo¬ 
graphic  interferometry).  On  the  other  hand,  composite  materials  can  form  special 
filters  for  use  in  acoustics. 

Biomechanics,  which  is  an  expanding  field,  can  supply  a  variety  of  natural, 
anisotropic,  composite  materials  (particularly  bones)  which  are  impervious  or 
porous.  In  return,  it  is  not  impossible,  in  the  near  future,  that  there  will  be 
a  development  in  the  use  of  composite  materials,  compatible  with  human  tissues, 
for  the  fabrication  of  prosthesis  or  orthesis  items. 

Finally,  it  should  be  noted  that  there  is  a  growing  im.  n 
matics  and  theoretical  mechanics,  in  composite  materials  which  are  vittiei 
periodic  or  which  have  random  distribution.  During  the  past  ten  years  there  has 
been  an  outburst  of  various  studies  on  such  media,  in  all  directions.  Mathe¬ 
matical  physics  methods  also  find  a  rich  field  there. 

Taking  into  account  the  multiple  aspects  which  composite  materials  offer, 
and  within  the  scope  of  the  present,  informal  report,  discussion  will  be  limited 
to  some  theoretical  and  experimental  studies. 

For  the  former,  two  aspects  will  be  examined,  assessment  calculations  and 
the  propagation  of  waves  and  vibrations  in  composite  materials.  This  done,  I  am 
very  conscious  of  leaving  aside  many  sectors  whose  importances  are  at  least  equal 
to  those  of  the  chosen  areas. 

2  ASSESSMENT  CALCULATIONS 

The  purpose  of  the  assessment  calculations  is  to  supply  elastic  or  visco¬ 
elastic  constants  of  composite  materials.  In  the  plastic  region  (or  at  breakage) 
they  allow  calculation  of  the  characteristics  of  the  boundary  surface.  The  idea 
is  to  replace  the  real  material  by  an  equivalent,  homogeneous  material. 

2 . 1  Boundary  methods 

Boundary  methods  have  given  a  variety  of  formulae  ranging  from  the  simplest 
to  the  most  complicated.  They  are  derived  from  known  theorems  relating  to 
elastic  potential  energy.  Two  approaches  are  used,  one  kinematically  allowable, 
the  other  one  statically  allowable.  If  the  geometry  of  the  medium  can  be  sacri¬ 
ficed,  then  the  simplest  formulae  are  those  of  Reuss  and  Voigt.  However,  when 
the  mechanical  characteristics  of  the  two  constituents  (or  phases)  of  the 
material  are  very  different  one  from  another  (elastic  constants  and  specific 
gravities)  the  boundaries  are  very  elongated  and  divergent  and  of  little  usefulness. 
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2.1.1  The  Hill  theory 

8  .  . 

leory  ,  arising 

applies,  particularly,  to  isotropic  materials  with  spherical 


8  9 

This  theory  ,  arising  out  of  the  work  of  Kroner  gives 


an  improvement, 
inclusions . 


It 


2.1.2  The  Hashin  and  Rosen  variational  methods 

These  authors*^  have  made  an  important  progress  in  the  field  of  assessment 
calculations,  concerning  materials  having  hexagonal  or  random  distributions  of 
unidirectional  fibres.  These  fibres  can  be  hollow  or  solid  (Fig  2). 


The  representative  elementary  volume  (,R-E.V.)of  the  orthotropic  material  is 
defined.  In  the  case  of  a  hexagonal  distribution,  the  admissible  fields  are 
those  which  produce  either  a  constant  deformation  or  a  constant  constraint  in  the 
space  contained  between  the  hexagon  (R.E.V.)  and  the  inscribed  circle  of  the 
cross  section.  The  problem  is  restored  to  a  classical  elasticity  problem,  in  a 
cylinder  of  radius  r.  with  limiting  conditions  of  the  type: 


u.  =£..x.  or  F.  =o..n. 

i  ij  J  i  1J  J 


Displacement  on  Surface  density 

the  boundary  S  of  forces 


(I) 


In  this  way,  Hashin  supplied  the  boundaries  for  the  five  elastic  constants 
in  the  case  of  hexagonal  symmetry  and  six  in  the  case  of  random  distribution. 

2.1.3  Variational  methods  using  polarisation  tensors 

Hashin  and  Shtri.kman!  *  have  made  a  further  improvement.  Instead  of  mini¬ 
mising  (or  maximising)  the  potential  energy  of  the  admissible  fields,  they 
minimise  (or  maximise)  a  certain  function  U  ,  defined  on  a  functional  space 
E(V)  ,  where  V  is  the  volume  and  E(V)  is  chosen  as  the  product  of  the  space 
of  kinematically  admissible  fields  and  the  space  of  polarisation  constraints  p.. 


o .  . 

ij 


Ci j  k£Ck£  +  Pij 


(2) 


where  a.,  and  c,  are  the  constraints  and  the  deformations  and  the  exponent  0 
ij  R2 

refers  to  the  comparison  medium. 

The  originality  of  the  method  lies  in  the  fact  that  U  is  not  defined 
uniquely  by  the  material  M  being  studied,  but  also  by  the  reference  material 
Mq  ,  having  the  same  geometry  and  being  subjected  to  the  same  boundary  conditions 
The  step  is  similar  to  that  adopted  in  electricity  for  the  study  of  dielectrics. 
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The  statically  admissible  fields  can  also  be  used,  together  with  a  polarisa¬ 
tion  displacement,  by  prescribing  a  force  field  on  S  . 

1  2 

Alblas  and  Kuipers  have  generalised  the  results  of  Hashm  and  Strikman  in 
a  space  Ej(v)  containing  E(v)  by  modifying  the  function  U  . 

Side  by  side  with  these  variational  formulae,  and  for  industrial  applica¬ 
tions,  there  are  many  assessment  methods,  some  semi-empirical  (Chamis  and 
.  1  3 

Sendeckyj  )  and  others  which  take  account  of  faults  in  the  composite  material 

1 4 

(fibre  contiguity,  imperfect  alignment,  etc),  such  as  that  by  Tsai 

2.1.4  Elastostatic  methods 

Chevalier*  \  dealing  with  the  hexagonal  fibre  network,  used  a  representa¬ 
tive,  elementary  volume,  as  did  Hashin.  He  proposed  the  replacement  of  the 
hexagonal  (R.E.V.)  contour  by  a  circular  one  which  maintains  the  percentage  of 
the  fibres  by  volume.  He  w^-ked  with  materials  having  transverse  isotropy. 

2.1.5  The  elastodynamic  method 

This  method  is  advocated  by  Behrens *^  as  being  useful  for  long  waves 
applied  to  a  variety  of  periodic,  composite  materials  (lamellar  composite 
materials,  composite  materials  with  anisotropic  fibres  or  with  sloping  elastic 
constants,  composite  materials  with  peculiar  symmetry,  etc).  Behrens  uses  plane 
waves  in  the  periodic  medium. 

The  methods  given  above  are  summarised,  for  comparison,  in  Table  I,  which 
is  taken  from  Chevalier'^.  Fig  3  and  Table  2  give  some  indications  on  three 
assessment  methods  applied  to  glass  fibre  composite  materials. 

To  conclude  this  section,  it  can  be  said  that  without  prejudice  of  knowing 
the  elastic  constants  of  each  phase  of  the  material  and  the  percentage  of  fibres 
by  volume,  there  are  now  available,  for  materials  with  transverse  isotropy, 
methods  of  assessment  calculation  which  a-e  making  their  mark  in  industrial 
applications.  In  my  opinion,  the  finite  elements,  in  spite  of  many  attempts, 
have  not  produced  any  significant  contributions  in  this  field. 

2. 2  Assessment  calculations  in  linear  viscoelasticity 

The  use  of  the  correspondence  principle  allows  the  viscoelastic  constants 
of  a  material  to  be  obtained  from  elasticity  formulae.  However,  existing  calcula¬ 
tions  are  only  true  in  the  context  of  perfect  adhesion  between  the  fibres  and  the 

matrix  and  do  not  take  into  account  the  diffraction  of  waves  at  the  interfaces. 

18  19 

On  this  subject  the  work  of  Hashin  and  Chevalier  should  be  noted. 
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2 . 3  Assessment  calculations  in  the  plastic  state  or  at  breakage 

.20  ...  .  . 

Kelly  and  Davis  have  used  them  for  unidirectional,  composite  materials. 

They  examined  the  following  three  distinct  processes: 

(a)  breakage  controlled  by  the  resistance  to  simple  fracture  along  the  fibre 
axes; 

(b)  transverse  resistance  in  a  direction  perpendicular  to  (a); 

(c)  pure  shear  resistance  of  the  matrix. 

2.3.1  The  MacLaughlin  approach 

Limit  analysis  has  been  used  with  more  or  less  success  by  various  research- 

21 

ers  into  multilayer  composite  materials.  MacLaughlin  and  Batterman  have  consid¬ 
ered  the  extreme  simplified  cases  (high  or  low  percentages  of  fibres)  which, 

unfortunately,  do  not  fall  within  the  normal  percentage  range  (60%  fibres). 

22 

Majumdar  and  MacLaughlin  have  attempted  to  state  the  general  principles  of  the 
application  of  limit  analysis  methods  in  the  search  for  a  boundary  surface.  Their 
equations  have  some  weaknesses. 

2.3.2  The  Le  Nizerhy  approach 

23 

In  particular,  the  work  of  Le  Nizerhy  should  be  mentioned.  He  has 
advocated  a  consistent  method  for  enclosing  the  site  of  the  mean  constraint 
vector  by  means  of  a  lower  envelope  (approach  through  constraints)  and  an  upper 
envelope  (dual  approach).  In  the  case  of  a  per f ec t 1 y-pl as t i c  rigid  system,  it 
requires  static  and  kinematic  approaches  to  the  boundary  surface  (classical  limit 
analysis).  Taking  into  account  the  fragile,  ductile  rupture  of  fhe  fibres, 

Le  Nizerhy  was  led  to  consider  the  calculation,  using  only  the  hypothesis  of 
convexity  as  criterion,  and  with  the  mathematical  formalism  unchanged. 

The  material  layout  takes  into  account  the  periodic  distribution  of  fibres 
both  in  each  layer  and  in  the  resin  which  forms  the  homogeneous  layers.  The 
loading  parameters  are  introduced  as  being  the  mean  constraints  (with  close 

homothetics)  in  a  macroscopical ly  homogeneous  volume  of  the  composite  material. 
The  deformation  velocity  vector  q  associated  with  the  vector  Q  in  the 
expression  for  the  principle  of  virtual  powers,  is  expressed,  this  latter  taking 
the  form: 


tr(o  .  d)d.U 

.  “ 


I 

*  1 


Qiqi 


=  V  (Si)  o  .  .d.  . 

i  ]  i.l 


(3) 
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From  Che  definitions  of  the  mean  quantities,  o..  and  d..  ,  can  be  obtained  an 
interior  envelope  with  potentially  supportable  loading  surfaces,  by  trying  to 
obtain  the  colvex  envelope  of  all  the  calculated  loadings  and  an  exterior 

envelope,  using  a  classical  method  of  approach  for  the  envelope  of  a  family  of 
planes.  Fig  4  shows  an  outline  of  how  the  method  can  be  used  in  an  application 
affecting  multilayer  composite  materials  with  boron  fibres. 


J  PROPAGATION  OF  WAVES  AND  VIBRATIONS  IN  COMPOSITE  MATERIALS 

It  is  in  this  field  that  a  large  number  of  theoretical  studies  are  to  be 
found.  In  order  to  grade  them,  the  first  thing  is  to  compare  the  wavelengths  A 
with  respect  tc  the  dimensions  of  the  elementary  network,  d  ,  of  the  composite 
material.  If  A/d  is  large,  the  behaviour  of  the  composite  material  is  that  of 
an  equivalent,  homogeneous  material  having  the  elastic  characteristics  given  by 
the  assessment  calculations.  In  this  context,  there  is  available  what  is  known 
as  the  theory  of  effective  moduli.  It  is  then  possible  to  tackle  the  study  of 
the  vibration  of  rods,  plates  and  shells,  on  the  one  hand,  and  the  propagation  of 
long  waves  in  composite  materials  on  the  other.  This  theory  finds  its  limitations 
very  rapidly  when  A/d  approaches  unity.  It  is  then  necessary  to  use  more 
elaborate  theories  which  could  take  better  into  account  the  dispersion  of  waves  at 
higher  frequencies.  In  this  context,  the  theory  of  continuous  media  with  micro¬ 
structure  appears  to  be  the  most  promising.  Another  approach  is  that  of  the  so- 
called  theory  of  mixtures  which  uses  the  interaction  between  the  components  of  the 
composite  material.  Variational  methods  which  affect  periodic  networks  have  been 

proposed  by  various  authors  and  have  given  rise  to  recent  interesting  improvements 
24 

by  Nemat-Nasser  .  The  theory  of  continuous  media  based  on  asymptotic  develop¬ 
ments  should  also  be  noted. 


3 .  i  Theory  of  effective  moduli 


3.1.!  Wave  propagation 


The  mean  constraints  a..  and  the  mean  deformations  e.. 

ij  _  iJ 

live  elementary  volume  are  defined  and  the  elastic  constants  of 


anisotropic,  homogeneous  medium,  are  calculated 


in  a  representa¬ 
tive  equivalent, 


o .  . 

i.l 


C*  r 
ijki  k£  ' 


(4) 


In  an  indefinite  medium,  the  study  of  wave  propagation  using  equation  (4) 
does  not  allow  the  <li s torsion  of  the  waves  to  be  taken  into  account.  However, 


4  J 


*?wz  n 
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this  theory  does  provide  a  ra-'id,  mechanical  characterisation  oi  the  medium  by 
ultrasonic  methods  using  long  waves. 

3.1.2  Vibrations  of  rods  and  plates 

A  certain  number  of  classical  vibration  studies  for  isotropic  rods  and 

plates  can  be  transposed  to  composite  materials.  For  rods,  these  include  rhe 

25  26  27  .28 

work  of  McNiven  ,  Le  Nizerhy-Vinh-Cheval ier  ,  Vinh-Nugues  and  Touratier 

relating  to  approximate  theories  on  rods. 

It  is  a  matter  of  extending  the  known  approximate  theories  of  Hermann- 

Mindlin,  of  Volterra  (longitudinal  vibrations),  of  Timoshenko's  method  (flexion 

29 

vibrations)  and  of  Barr-Engstrom' s  method  (torsional  vibrations). 

It  should  be  noted,  in  passing,  that  in  general  the  dispersion  of  waves  for 
the  three  modes  of  decoupled  vibrations  is  much  more  pronounced  for  composite 
materials  than  for  isotropic  ones. 

This  dispersion  is  due  to  the  strongly  anisotropic  character  of  the  medium 
studied,  which  is  shown,  par ticular ly ,  in  the  high  ratio  between  the  longitudinal 
Young's  modulus,  E^  ,  and  the  shear  modulus  p  : 

Ej/p  as  20  to  50,  instead  of  about  3,  for  isotropic  materials. 


Here  and  now,  for  the  choice  of  displacement  fields,  polynomials  of  a  high 

degree  must  be  used  to  translate,  correctly,  the  dispersion  of  phase  and  group 

27 

velocities.  Applying  dynamic  torsion  to  beams  of  rectangular  cross  section*"  , 
the  Saint-Venant  theory,  which  relates  to  the  buckling  function,  reaches  its 
limits  in  the  study  of  the  higher  modes. 

Rods  (whose  axes  are  not  confounded  with  axes  of  symmetry  of  the  material) 

give  rise  to  coupled  waves.  Several  approximate  theories  exist,  such  as  those  by 
30  .  31 

Abarcar  and  Cunmff  and  Le  Nizerhy  .  They  all  require  complementary  studies. 

.32 

The  method  advocated  by  Mindlrn  of  the  double  series  for  the  displacement 
’omponents  should  be  noted. 


II 


m  n  m, n  . 
X2X3Ui  (X1 


t) 


I*' 


m  n  m,  n 
2X3Ui 


(5) 


where  x^  is  the  beam  axis  and  t  is  time, 
tion  (51  depends  upon  the  modes  of  vibration 


The  number  of  terms  taken  in  equa- 
studied.  For  plates,  the  well-known 


I  2 


approximate  theory  by  Mindlin  has  been  extended  to  composite  materials  by 
33 

Whitney 

Amongst  the  work  in  France  should  be  noted  that  of  Verchery^  on  composite 
plates  in  a  static  state.  This  work  is  capable  of  being  extended  to  the  dynamic 
state. 


3  .  - 


Theories  of  continuous  media  with  discrete  structures 


A  number  of  theories  with  microstructures  have  been  proposed,  amongst  which 

34  33  36 

are  me  works  of  Mindlin  and  Tiersten  ,  the  theories  of  Toupin  and  Truesdell 

37 

and  the  so-called  micropolar  theory  of  Er ingen  .  This  latter  tends  to  general- 

38 

ise  the  other  theories.  A  comparative  study  by  Eringen  brings  out  the 
differences  between  these  theories  (restrictions  imposed  on  the  spatial  micro¬ 
rotation  vector  ±,  =  1< 


"kim  m, H 


where  u  is  the  displacement  vector;  the 
m 


presence  of  the  couple  per  unit  mass  and  of  the  acceleration  u^  in  the 

behaviour  laws;  omission  of  the  micropolar  rotatory  inertia  term;  final  depen¬ 
dence  of  the  constraint  vis-^-vis  the  applied  loads  and  the  inertia,  etc). 

Most  of  the  theories  on  microstructures  are  concerned  with  linear, 
isotropic  materials. 

As  regards  the  propagation  of  waves  in  an  indefinite  medium,  the  following 
theoretical  results  by  Eringen  should  be  noted. 

For  plane  waves,  Fig  5  shows  that  in  the  direction  of  propagation,  n  , 
there  is  a  microrotational  plane  wave  in  addition  to  the  classical  longitudinal 
wave.  For  the  microrotational  wave,  propagation  is  only  possible  if  the 
frequency  exceeds  a  certain  cut-off  frequency.  Fig  6. 


For  transverse  waves,  in  addition  to  the  classical  transverse  displacement 

wave,  U  ,  Eringen  has  found  a  wave  with  microrotation  if  which  is  orthogonal  to 

both  th .  direction  of  propagation,  n  ,  and  the  displacement  U  .  This  latter 

transverse  wave  with  microrotation  only  exists  for  frequencies  greater  than  a 

critical  frequency.  Below  this,  it  degenerates  into  an  evanescent  wave  attenua- 

38 

ted  by  distance.  For  surface  waves,  Eringen  and  Suhubi  have  predicted  a  second 

type  of  dispersive  wave.  Here,  too,  no  experimental  confirmation  has  yet  been 

put  forward.  Amongst  the  theories  for  discrete  structures  applicable  to 

composite  materials  should  be  mentioned  two  proposals,  one  by  Hermann,  Sun  and 
39 

Achenbach  ,  known  as  the  effective  stiffness  theory  and  the  other,  the 
4  I 

Drumhel ler-Bedford  theory  which  uses  a  second  order  approximation. 
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3.2.1  Theory  of  effective  stiffness 

This  was  developed,  originally,  to  interpret  dynamic  effects  such  as  the 
dispersion  of  waves  in  stratified  media  (multilayers).  Fig  7  represents  a 
rectangular  disposition  of  unidirectional  fibres  in  a  plane  cross  section.  In  a 
grid  (k,  Z)  it  is  a  matter  of  adopting  series  for  displacements.  The  following 
developments  can  be  adopted  fairly  well  in  the  fibre  (f)  and  matrix  (m) : 


f(k,0 

-(k,A) 

.  f (k,  £)  x  .  . 

. f (k, Z) 

i . 

i 

=  u . 
i 

+  r  cos  0  .  ip..  +  r  sm  0  , 

2i 

'  ^3i 

m(k,  t) 

-(k,  £) 

f (k, £) 

m(k, £) 

i. 

l 

=  u . 

l 

+  a  cos  0  .  ib..  +  a  sm  0  , 

2i 

■hi 

(fibre) 


/  %  „  ,m(k,£)  .  .  ,m(k,£)  ... 

+  (r  -  a)  cos  0  .  i +  (r  -  a)  sm  0  .  (6) 


with  r  >  a  matrix  . 


-f(k  Z) 

The  physical  interpretation  is  simple,  u.  ’  represent  overall  displace- 
m(k  1)  ^  . 

ments  and  ’  (n  =  2,  3)  represent  local  displacements.  The  two  types  of 

displacements  are  defined  for  discrete  values  of  x 2  and  and  are  continuous 

functions  of  the  variable  x^  and  time.  Since  it  is  not  possible  to  satisfy, 
exactly,  the  conditions  for  continuity  of  the  displacements  at  the  interfaces 
between  neighbouring  cells.  Sun  and  Achenbach  suggest  a  continuity  for  the  mean 
displacements. 

On  the  basis  of  these  series  developments  (which  use  Legendre  polynomials) 
they  have  evaluated  the  kinetic  and  potential  energies  whose  means  are  defined  at 
the  fibre  centres.  The  transition  from  a  discrete  lattice  to  a  continuous, 
homogeneous  model  is  accomplished  by  defining  the  fields  for  the  kinematic  and 
dynamic  variables  which  are  continuous  for  all  the  coordinates.  Thus  the  energy 
deformation  density  contains  not  only  the  effective  moduli  terms,  shown 
previously,  but  also  terms  with  constants  arising  from  combinations  of  the 
geometric  and  elastic  properties  of  each  phase.  The  equations  are  obtained  by 
the  Hamilton  principle  used  in  conjunction  with  the  continuity  relationships  and 
Lagrange  multipliers. 

As  a  particular  case,  the  researchers  have  studied  multilayer  plates  where 
each  layer  is  isotropic.  Fig  8  shows  the  effective  stiffness  theory  with  a  first 
order  approximation.  It  shows  that  for  the  phase  velocity  agreement  is  good  for 

very  low  modes  (number  of  wave  k  <  2  with  respect  to  the  exact  theory  of 

„  40. 

Rytov  ) . 
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The  model  with  a  first  order  approximation  translates,  correctly,  the 
dispersion  of  shear  modes  but  not  the  other  modes.  A  second  order  approximation. 
Fig  9a,  gives  a  better  approximation  to  the  exact  theory. 

3.2.2  General  theory  of  effective  stiffness 

As  an  application  of  the  previous  theories,  the  problem  of  a  semi-indefinite 

or  finite  body  subjected  to  a  permanent  or  transitory  displacement  and/or  the 

limiting  constraint  conditions,  is  of  ^reat  importance.  These  problems  have  been 

41 

formulated  by  Bedford  and  Drumheller  They  can  be  applied  to  the  highest  order 


of  effective  stiffness  theories.  These  two  authors  introduce  an  external,  plane 
surface  and  systems  of  coordinates  x^,  (Fig  9b)  relative  to  a  multilayer 

composite  material,  as  well  as  local  coordinates  (x^,  xM  and  (n,  s)  whose 
origins  are  situated  on  the  fibre  axis.  The  displacement  vector,  ,  in  the 

kth  layer  can  then  be  written  as: 


fk  -fh/-f  _f  \ 
ui  =  Ui  \WC)  = 


«fh,_  _  , 

u.  (n, s , t) 


(7) 


The  limiting  conditions  for  the  displacements  in  terms  of  the  distribution  of 

fk 

specified  displacement  IK  (s,  t)  at  the  external  boundary  of  the  layer  k  are: 


-fk/_f  _f 
u.  lx.,x„ 
l  \  i  2 


’  C)  n=C 


*fk ,n  -  , 

(0, s, t) 


..fk,-  . 

U.  (s,t) 


(8) 


a  f  lc  f  k  .  _ 

Then  and  are  developed  in  terms  of  s  and  the  development  is 

of  the  same  order  as  the  development  of  displacement  vectors  in  the  theory 

adopted . 


By  equalising  terms  in  the  same  power  of  s  ,  a  set  of  limiting  conditions 
is  obtained: 
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This  procedure  is  direct  and  shows  the  relationship  between  the  limiting  condi¬ 
tions  and  the  distributions  of  displacements  and  constraints  over  the  external 
surface . 

The  problem  of  the  accuracy  of  the  effective  stiffness  theory,  with  its 

4  1 

limiting  conditions,  has  been  discussed  by  Bedford  and  Drumheller  .  They  have 
shown  that  the  microstructure  theory,  with  a  finite  number  of  degrees  of  freedom, 
has  only  limited  possibilities  for  modelling  the  distributions  of  displacements 
and  constraints  at  the  boundary. 

3 . 3  So-called  theory  of  mixtures,  or  of  media  with  interactions 

This  is  a  matter  of  a  theory  used,  originally,  in  geology  and  recently 
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adapted  to  composite  materials  by  Lempri&re  ,  Bedford  and  Stern  and  Hegemier 

It  is  assumed  that  each  phase  of  the  composite  material  is  subjected  to 
individual  deformations.  The  microstructure  of  the  composite  material  is  shown 
by  the  interactions  between  the  constituents.  Some  indications  of  how  this 
method  is  applied  to  multilayers  will  be  given  here. 

By  integrating  the  equations  of  movement  in  the  deformation  plane  we  have: 


ij  -J 


=  pu 


(a) 
i ,  £t 


(10) 


where  a  indicates  one  of  the  two  phases  of  the  medium,  with  respect  to  the 


variable  x 


(a) 


in  the  direction  of  the  thicknesses  h 


(a) 


obtained  of  the  form: 


Equations  are 


(a)  3°ll)e 


3x 


h(a)  (a) 
h  p 


32  (a) 

d  u,  e 


3t 


(a) 

13 


(x),h(a)t) 


(11) 


where  the  exponent  e  indicates  a  mean. 


r(a) e ' 


h(a)  /(a)' 


,(a) 


dx 


(a) 


(12) 


••  (&)•«  •• 
and  assures  the  continuity  of  the  term  t._  ,  which  is  asymmetrical  with  respect 

to  *$«>  : 
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(0  (1) 


(2)  (2) 


t  j  3 (x  j , h,  t)  =  T]3(x,h,t)  =  T*3(xt,t) 


Equations  of  movement  are  obtained,  of  the  form: 


(!)  (!)  (I) 


3o 


1  1 


3x, 


-  P 


3t2 


Pc3 


(13) 


(14) 


where  the  exponent  p  indicates  partial  constraints  or  displacements. 


(a)  (a)(a) 


D  C 

ail  =  V  GH 


(a)  (a)  (a) 


=  V  p 


(15) 


(a)  (a)(1)  (2) 

V  =  h/(h  +  h)  , 


(a) 

where  V  is  the  volume  fraction, 


(1)  (2) 

p  6  =  t*  /(h  +  h) 

c  I  J 


(16) 


with 


(a) 


(a) 

c 


(17) 


Thus  the  p  6  term  in  equation  (14)  represents  the  transfer  of  the  amount 
c 

of  movement  from  one  phase  to  another.  Hence  it  corresponds  to  a  volume  force 
equivalent  to  an  interaction  due  to  shear  at  the  interfaces. 

The  basic  problem  is  the  determination  of  these  interaction  terms  p  6  and 

c 

the  equations  of  behaviour  for  each  phase  of  the  composite  material.  For  this, 

44 

Hegemier  et  ai  have  suggested  an  asymptotic  development  as  a  function  of  x3 
for  the  constraints  and  displacements. 


The  equations  of  behaviour  are  of  the  form: 
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(18) 


which  express  the  constraints  in  terms  of  the  composite  material  constituents  and 
the  displacement  gradients  in  each  phase. 

Figs  10  and  11  show  the  phase  velocities  in  terms  of  the  frequency  for  two 
types  of  multilayer  composite  materials.  The  agreement  between  theory  and 
practice  is  good  for  the  low  frequencies. 

3 . 4  Theory  of  continuous  media  with  microstructure  based  on  asymptotic 
developments 

Several  formulae  exist  for  the  propagation  of  waves,  such  as  those  of 
47 

Hegemier  and  of  Nayfeh 

The  method  given  by  Nayfeh  is  of  particular  interest  when  the  number  of 
layers  exceeds  three  in  each  cell  and  the  pattern  is  repeated  particularly  in  one 
direction.  For  lack  of  space,  only  the  steps  adopted  for  a  unidirectional 
problem  will  be  given  here. 

(a)  As  in  the  theory  of  mixtures,  the  equation  of  movement  is  integrated  with 
respect  to  the  variable  x  in  each  elementary  cell  (0  <  a  <  m) . 

(b)  By  assuring  the  conditions  of  continuity  between  one  cell  and  another,  the 

(k)  (k) 

relationships  are  sought  which  give  the  constraints  o^  (0)  and  u^  (0)  ,  the 
displacements  at  the  median  plane  of  each  layer. 

(c)  Next,  the  constraints  and  displacements  found  above  are  developed  as 

•  •  (k)  (k) 

series,  as  functions  of  the  local  coordinate  y  =  x  -  x  ,  where  x  is  the 

(k)  “  .  a 

overall  coordinate  and  x  is  the  local  coordinate  in  the  a  layer  of  the  kth 

a 

cell.  Thus  a  range  of  finite  difference  equations  is  obtained  for  the  spatial 
and  differential  coordinates  with  time. 

(d)  Then  a  small  parameter  e  =  A /£,  is  introduced,  where  A  is  the  total 
thickness  and  l  is  a  macrodimension  of  the  problem,  which  could,  for  example, 
represent  the  wave  length  A  .  When  e  -+  0  the  phase  velocity  obtained  is  that 
given  by  the  theory  of  homogenised,  non-dispersive  media. 
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Amongst  the  approximate  theories  put  forward  up  to  now  with  models  of  wave 
dispersion  for  the  lowest  modes,  the  method  of  asymptotic  development  gives  a 
greater  accuracy  than  do  the  others. 

3 . 5  Variational  methods 
3 

Kohn  has  applied  the  theory  of  Floquet,  Bloch  and  Bnllouin  to  composite 

48  ...  .  . 

materials.  Lee  has  recently  reviewed  various  variational  methods  in  detail. 

The  variational  principles  are  developed  in  an  integral  form  in  a  cell  which 

represents  the  representative  elementary  volume.  The  deformation  energy  and  the 

complementary  energy  have  been  used.  In  passing  from  the  fibre  to  the  matrix, 

the  discontinuities  in  certain  constraints  and  deformations  at  the  interfaces  are 

taken  into  account.  A  Rayleigh-Ritz  procedure  has  been  applied  to  calculate  the 

dispersion  relationships  for  any  propagation  direction  whatsoever.  So  far  as  is 

known,  the  studies  by  Lee  have  been  confined  to  multilayers  (periodicity  in  one 

direction).  The  credit  for  a  generalisation  of  such  methods  and  for  an  applica- 

24  49 

tion  to  several  cases  rests  with  Nemat-Nasser  ’  who  has,  in  addition,  proposed 
a  dual  Rayleigh  coefficient  (based  on  the  constraints).  For  various  situations 
(pronounced  density  discontinuity,  discontinuity  of  the  elastic  constants)  these 
variations  of  the  coefficients  show  the  Rayleigh  coefficient  to  be  adequate.  A 
new  Rayleigh  coefficient,  based  on  mixed  fields,  appears  to  be  very  effective  and 
gives  the  upper  boundaries  for  particular  values  with  needing  recourse  to  the 
usual  or thogonalisation  of  the  test  functions. 

3.5.1  The  Rayleigh  displacement  coefficient 

This  coefficient  is: 


X 


R 


(  pu,u  > 


(  nu' ,u’ ) 
(  pu,u> 


1 

J'  pu  .  v*dx 

I 


> 


(19) 


where  p 
u 

n 

* 


is  the  density, 

is  the  displacement  in  reduced  variable, 

is  the  elastic  coefficient, 

indicates  the  conjugated  complex  quantity, 

indicates  differentiation  with  respect  to  x  . 
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The  periodicity  of  the  medium  is  expressed  as: 


u(|) 


u(-j)e 


iQ 


(20) 


where  Q  is  the  wave  number. 

If  L  is  positive,  auto-associated  operators,  for  harmonic  waves  we  have: 

Lu  =  -  -y.--y  [  n(x)u'  J  =  Au  -*  <  x  <  i  (21) 


At  the  points  x. 
following  conditions: 


(j  =  1,2) 


of  discontinuity  n  ,  we  should  have  the 


n(x.)u'(x.) 
J  J 


n(xt )u' (xT) 


x.  =  lim(x.  ±  e) 
J  e-K)  J 


j  =  1  or  2 

e  -*■  0 


(22) 


Equations  (20)  and  (21)  can  be  satisfied  by  taking  for  the  displacement: 


-M ' 

exp  i(Q  +  2iroc)x 

a=0 


(23) 


The  method  of  calculation 


is  the  classical  one: 


3X(5) 

3U* 

a 


0,±1,...,±M'  .  (24) 


When  n  has  discontinuities  or  varies  very  rapidly,  this  method  is 
ineffective . 

3.5.2  The  dual  Rayleigh  coefficient 
Assuming  that: 


we  have 


R 


l 

P 


and  Do  =  u' 


(25) 


(Ro')  +  ADo  =  0  , 


-J  <  x  <  i  (26) 


i 


AR  is  made  stationary: 


3Xr(5) 

3S* 

a 


=  0 


a  =  0,±1 . -M' 


When  R  =  1/p  is  not  continuous,  AR  is  not  effective. 
3.5.3  A  new  coefficient 

Taking  the  mixed  fields  a  and  u  we  obtain: 


<q,u')  +  (u',o)  +  (00,0) 

<  DU,V> 


(29) 


It  can  now  be  shown  that: 


and  that 


AR(u)  >  XN(u,o) 


=  A.^  >  A  ,  if  n  is  a  constant 

R  N  p 


(30) 


3.5.4  A  new  dual  coefficient 


wi  th 


and  if  p  is  a  constant: 


<V,o')  +  (q',W)  -  (pW,W) 
<  Do , o  > 


Ar(«j)  >  An(o,W) 


(31) 
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3.5.5  Improvements  in  the  new  coeff icients 

This  has  led  Nemat-Nasser  to  propose  a  modification  to  the  test  functions 
u  and  a  ,  as  in  equation  (22),  so  that  the  following  continuity  condition  is 
not  violated. 


R(x . )  .  a(x.) 
J  J 


R(x+)  .  o' (x+) 
J  J 


with  j  =  1  or  2  .  (32) 


Moreover,  he  ensures  that  the  new  test  functions  u  and  a  are  such  that 
qu '  and  R'V  are  continuously  differentiable  in  [—3,5  ]. 

For  the  unidimensional  problem  the  modified  functions  are: 


±M' 

s  '  I 

a=0 


(exp  iQ  —  1  ) 


“1 

2 


D(C)  exp  i  (Q  +  2t!a)  ddf  +  j  D(£)  exp  i  (Q  +  2iTa)  £d£ 

. (33) 


-M ' 

■5'. 


(exp  iQ  -  1 ) 


p(C)  exp  i(Q  +  2ta)  CdS  +  j  p(£)  exp  i  (Q  +  2trciKdf; 

“  2 

. (34) 


,  and  a  new  upper  boundary, 
R 


,  are  obtained  at 


and  a  new  lower  boundary, 
the  price  of  an  improvement  in  u  and  0  . 

Nemat-Nasser,  Fu  and  Minagawa^  have  also  provided  boundaries  for  the 
natural  frequencies  in  three-dimensional  problems.  These  researchers  have  made 
numerous  applications  of  the  method,  which  indicates  the  interest  shown  in  it. 

Mention  should  be  made  of  the  work  of  Gotteland^*  who  has  carried  out  a 
systematic  study  of  the  various  variational  formulae  with  a  view  to  the  applica¬ 
tions  to  multilayer  plates  by  finite  elements. 


3.6  Various  models  which  take  wave  characteristics  into  account 


The  propagation  of  waves  in  composite  materials  is  always  accompanied  by  a 
attenuation.  This  can  be  explained  by: 

(a)  the  periodic  character  of  the  lattice; 

(b)  the  viscoelastic  properties  of  the  constituents; 

(c)  the  diffraction  and  diffusion  of  waves  by  the  fibres,  by  interfacial 
defects,  by  holes  in  the  matrix,  by  bad  fibre  alignment,  etc. 


1 
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Taken  as  a  whole,  the  problem  is,  theoretically,  of  great  complexity.  In 
practice,  there  have  been  only  tentative,  partial  explanations  of  certain  aspects 
relating  to  wave  attenuation. 

3.6.1  Modelling  by  analogy  with  a  continuous  medium 
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This  has  been  suggested  by  Drutnheller  and  Sutherland  ,  using  the  results 
obtained  by  Brillouin'  for  discrete  lattices  based  on  the  Floquet-Bioch  theorem 
relating  to  the  propagation  of  a  harmonic  disturbance.  The  wave  has  a  periodic 
dispersion  spectrum.  The  frequency  curve,  in  terms  of  the  wave  number  R  =  2tt/A  , 
where  A  is  the  wavelength,  has  frequency  bands  in  which  no  wave  can  be  propaga¬ 
ted  and  for  a  given  value  of  f  there  is  an  infinite  number  of  values  of  k  . 

The  classical  dispersion  relationship  can  be  used  for  the  lattices,  giving: 


„  ,  ,  (Z  +  l)2  .  ,  „  +  ,  (Z  -  l)2 

cos  (kh)  =  - — -  cos  (wt  ) - 2^ - 


cos  (ojt  t  ) 


where  Z  is  the  ratio  of  the  acoustic  impedances  of  the  two  layers, 
t+  is  the  sum  of  transit  times  in  the  elementary  cell, 

t  is  the  ratio  of  the  transit  time  difference  to  the  sum  of  transit  time, 

r 

h  is  the  total  thickness. 

By  analogy  with  a  continuous  medium,  the  properties  of  this  latter  can  be  used, 
and  are  as  follows: 

(i)  frequencies  filtered  in  the  by-pass  filter; 

(ii)  at  almost  zero  frequencies,  the  lattice  and  the  composite  material  have  the 
same  longitudinal  phase  velocity.  Then  equation  (35)  is  developed  and 
gives  the  relationship  for  the  dispersion  of  waves  from  the  lowest  mode; 

(iii)  a  redistribution  of  masses  is  undertaken  in  the  lattice  to  adapt  the 
dispersive  characteristics  of  the  ideal,  periodic  lattice  to  the  composite 
material . 

This  method,  which  is  very  simple  to  apply,  can  be  used  to  study  shock 


3.6.2  The  Sve  method 

The  Fourier  integral  has  been  used  to  simulate  the  remote  field  of 
constraints  near  to  shock  wave  fronts  by  assuming  that  the  predominant  modes  are 
those  with  the  lowest  frequencies.  Dispersion  is  included  by  adopting  a  complex 
wave  number.  Sve  studies  the  differences  and  the  common  points  between  the 
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geometric  dispersion  and  the  spatial  attenuation  of  the  waves.  He  uses  the  Airy 
function. 

3.6.3  Interpretation  of  dispersive  effects  by  the  Boltzmann  relationship 
54 

Christensen  has  suggested  the  use  of  dielectric  theory  as  a  possible 
approach,  since  the  two  basic  mechanisms  of  dielectrics  have  a  certain  analogy 
with  those  of  composite  materials: 

(i)  there  is  the  loss  due  to  displacement  of  loads  which  are  connected, 
elastically,  to  an  equilibrium  position; 

(ii)  there  is  the  loss  due  to  the  transition  of  loads  or  dipoles  between  the 
equilibrium  positions,  separated  by  a  potential  barrier. 

According  to  Christensen,  case  (i)  can  give  a  direct  character istation  of 
elastic  media  with  periodic  layers  and  case  (ii)  has  been  transposed  to  a 
homogeneous,  viscoelastic  medium. 

In  this  way,  Christensen  has  used  this  analogy  to  re-examine  the  dispersion 
relationship  obtained  by  Kohn  in  the  form: 

2  2.  2  4 

w  =  CqIc  -  6k  (36) 

where  and  6  are  parameters  which  express  the  layers  and  their  geometries 

in  terms  of  the  elastic  characteristics. 

He  has  also  examined  the  random  medium  by  adopting  a  relaxation  function  of 
the  form: 

G(t)  =  Ge"t/T  +  GQ 

where  G,  t  and  G  are  unknown  parameters  which  are  to  be  determined.  He 
obtains  the  attenuation  in  the  form: 

Gt)  CoWC,.Gl  '  P2G2)2 

n  -  - i - -O -  <37> 


where  C 


max 

V  P2 

G,,  G2 
h1  ’  h2 


is  the  maximum  wave  velocity, 
are  densities, 

are  the  moduli  of  elasticity  of  the  layers, 
are  the  thickness, 
is  the  pulse  rate. 
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Equation  (37)  allows  the  order  of  magnitudes  of  the  attenuations  of  waves  in 
magnesium  and  aluminium  to  be  remeasured. 

3.6.4  Diffraction  and  diffusion  of  elastic  waves 

As  far  as  I  am  aware,  very  little  study  has  been  given  to  this  problem, 

when  its  complexity  is  taken  into  account.  It  is  a  matter  of  studying  the 

diffusion  of  waves  due  to  the  frequency  of  fibres  sunk  in  a  matrix.  Bose  and 
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Mai  have  examined  the  case  of  the  transverse  wave  being  propagated  at  right 
angles  to  the  fibres  and  having  the  parallel  polarisation  of  the  fibres.  It  is  a 
case  of  random  fibre  distribution. 

These  authors  have  been  able  to  evaluate  the  complex  wave  number  of  the 
composite  material.  The  imaginary  part  of  the  wave  number  is  related  to  the 
geometry  of  the  fibre  lattice.  They  have  been  able  to  calculate  the  specific 
absorption  in  terms  of  the  fibre  concentration. 

4  SOME  EXPERIMENTAL  STUDIES 

4 . 1  Measurements  of  elastic  constants  in  micropolar  elasticity 
4.1.1  Static  measurements 

It  has  been  shown,  earlier,  that  certain  microstructure  theories  have  been 
used,  successfully,  in  the  study  of  wive  dispersion. 

A  serious  problem  arises,  namely  how,  within  the  scope  of  these  theories, 
can  the  elastic  constants  of  materials  be  evaluated  in  sufficient  numbers? 
Experimentally,  this  question  underlines  some  difficulties.  In  research  into  the 
laws  of  behaviour  in  micropolar  elasticity,  certain  attempts  have  been  made  which 
merit  further  notice. 

Thus  Gauthier  and  .lashman^  have  made  measurements  on  an  artificial,  micro- 
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lar  medium  tunned  1 

com  .i  1 

umi nium  particles  in  an  epoxv  resin  matrix. 
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where  a..  is  the  constraint, 
ij 

e..  is  the  deformation, 
ij 

i .  is  the  microrotation  vector. 
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m.  .  is 
ij 

X  and 
X,  a,  6 


the  restraining  couple, 

p  are  the  classical  Lam£  moduli,  and 

and  y  are  the  form  constants  inherent  in  the  micropolar  theory 
to  be  determined. 


For  an  isotropic  body  six  elastic  constants  have  to  be  evaluated.  From  equations 
(38)  and  (39),  X  has  the  dimensions  of  pressure  whilst  a,  6  and  y  have  the 
dimensions  of  forces  (couple/length). 

The  measurements  can  use  well-known  methods  of  classical  elasticity 
(traction,  torsion,  shearing,  circular  flexion,  etc)  for  rods  and  plates.  The 
analysis  formulae  are  more  complex,  for  example  in  the  case  of  traction  which 
gives  Young's  modulus  and  the  Poisson  number  in  the  following  forms: 


E  =  (2p  +  X) (3A  +  2p  +  X) (2X  +  2u  +  X) 


(2X  +  2p  +  X) 


(40) 


For  torsion  and  flexion  it  has  been  possible  to  demonstrate  the  character¬ 
istic  lengths  which  allow  a  definition  of  the  limits  within  which  micropolar 
elasticity  measurements  hold  good. 

Gauthier  has  shown  that  in  resolving  three  problems  at  the  limiting 
conditions  he  has  been  able  to  evaluate  the  six  constants. 

For  shear  measurement,  use  can  be  made  of  holographic  interferometry.  For 
torsion,  the  displacements  can  be  measured  by  double  exposure  hologram  and  a 
strongly-coherent  laser  beam.  These  optical  methods  avoid  physical  contact  and 
hence  local  disturbances  in  the  areas  of  interest. 

4.1.2  Dynamic  measurements 

Certain  methods  of  physical  measurement  can  be  used  to  evaluate  the  elastic 
constants.  A  list  of  these  methods  is  shown  in  Table  3. 

In  general,  the  optical  methods  given  are  concerned  with  short  wavelengths 
and  can  be  used  when  the  characteristic  micropolar  wavelengths  are  of  the  same 
order  of  magnitude  as  the  dimensions  of  the  lattice  (=  10  ^  cm).  It  should  be 
noted  that  Raman  spectroscopy  is  reserved  for  crystals. 

As  Askar"^  has  so  rightly  noted,  there  is  a  large  quantity  of  experimental 
results  on  crystals  which  await  interpretation  by  theor iticians  and  experimental 
workers  in  the  field  of  continuous  media. 
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Amongst  the  experimental  methods  which  can  show  up  certain  'micropolar 
effects'  (Table  2),  ultrasonics  and  hypersonics  are  very  serious  candidates. 

4 . 2  Mechanical  characteristics  of  anisotropic  composite  materials 

The  techniques  for  measuring  the  elastic  constants  of  non-polar  elastic  or 
viscoelastic  bodies  tend  to  increase.  Particularly  in  the  aeronautical  field 
there  are  important  developments  in  static  or  dynamic  experimental  techniques. 
Rods,  plates  and  tubes  have  been  used  to  this  end.  It  should  be  mentioned  that 
the  measurement  of  certain  'non-diagonal'  elastic  constants  (that  is  to  say  which 
are  not  situated  on  the  diagonal  of  the  corresponding  matrix)  avoids  couplings 
(f lexion-tor sion,  traction-torsion,  etc)  in  quasi-static  tests.  If  ultrasonics 
are  used,  it  will  be  possible  to  use  coupled  waves  whose  use  is  more  delicate 
than  decoupled  waves,  as  Curutchary'^ ,  Vinh  and  Garceau"^  and  Markham*^  have 
shown. 


Fig  12  shows  an  elasticimeter  with  alternate  flexion  and  torsion,  which 

allows  Young's  moduli  (E^),  shear  moduli  (p„)  and  the  coupling  constants  of 

composite  matei  ials  to  be  determined  on  samples  of  rectangular  (or  circular) 

cross  section.  Taking  into  account  the  influence  (which  is  not  negligible)  of 

the  rotational  inertia  effect  (in  flexion)  and  the  shearing  effect,  the  equations 

2  6 

of  motion  are  of  the  Timoshenko  type  with,  possibly,  coupling  terms  .  In  alter- 
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nate  torsion,  for  rectangular  cross  sections,  warping  must  be  taken  into 
account . 

In  other  respects,  the  use  of  samples  of  finite  length  complicates 
research  into  the  natural  modes,  even  by  small  amounts. 

58 

Fig  13  shows  some  experimental  equipment  which  uses  ultrasonics  .  Fig  14 
shows  the  slowness  (inverses  of  phase  velocities)  diagrams  relating  to  a  uni¬ 
directional  Kevlar  (fibre  in  aromatic  polyamide)  composite  material.  It  can  be 
seen  that  for  an  angle  of  incidence  of  more  than  13°,  the  longitudinal  wave  is  no 
longe’"  transmitted.  The  composite  material  therefore  plays  the  role  of  a  very 
directional  spatial  filter.  The  use  of  coupled  waves  for  the  evaluation  of  the 
last  'non-diagonal'  constant  gives  rise  to  Fig  15.  The  error  perpetrated  on  the 
value  of  the  constant  can  have  a  strong  influence  on  the  value  of  Cji^  • 

4 . 3  Study  of  wave  dispersion  by  ultrasonics 

Propagation  of  expanding  ultrasonic  waves  perpendicular  to  the  layers  has 
shown  both  waves  of  low  mode  (acoustic  branch)  and  those  of  a  higher  mode  (so- 
called  optical  branch).  Fig  16  shows  that  the  composite  material  behaves  like  a 
cut-off  filter  in  certain  frequency  intervals.  The  effect  of  this  property  is  to 


LT  2044 


LT  2044 


27 


r- 

;  I 

< 


modify  the  frequency  spectrum  of  the  progressive  wave,  as  is  shown  in  Fig  17 
(from  Robinson  and  Lippelmeier^ * ) . 

The  wave  dispersion  can  be  shown,  easily,  by  using  several  transducers 
having  a  range  of  natural  frequencies. 

4 . 4  Degeneration  of  ultrasonic  waves 

A  decoupled  wave,  applied  to  a  composite  material,  and  being  propagated 

along  a  direction  of  symmetry  of  the  material,  can  be  transformed  into  coupled 

waves  (transversal  and  longitudinal  waves).  This  happens  if  the  fibres  undergo 

directional  variations.  Hence  this  wave  transformation  can  be  used  to  study 

6  2 

certain  faults  in  fibre  alignment  (Weistman  ). 

4 . 5  Geometric  and  viscoelastic  dispersion 

Even  in  an  indefinite  medium  formed  from  unidirectional  fibres,  the 
geometrical  dispersion  of  the  waves  can  be  explained  by  fibre  size  (which  behave 
as  wave  guides)  and  by  their  periodicity.  The  viscoelastic  type  of  dispersion  is 
due  to  the  matrix  (often  of  resin).  These  two  types  of  dispersion  exist  together 

and  can  either  act  in  the  same  sense  or  can  neutralise  each  other,  depending  upon 

,  58,59 

the  modes 

4 . 6  Elastic  waves  in  limited  media 
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The  study  of  surface  waves  ,  of  plate  waves,  of  Love  waves  (surface  waves 
with  transverse  polarisation  which  are  propagated  between  a  substrate  and  a  layer) 
and  of  Stoneley  waves  (interface  waves)  has  not  given  rise,  as  far  as  I  am  aware, 
to  theoretical  and  experimental  studies  on  composite  materials.  The  studies 
would  certainly  be  more  refined  than  those  done  on  homogeneous,  isotropic  media 
but  they  could  lead  to  very  interesting  applications,  particularly  in  the  study 
of  adhesion  between  the  layers  of  composite  materials. 

4 . 7  Evidence  of  special  waves  related  to  the  micromorphic  nature  of  the  medium 

Section  3.3  drew  attention  to  particular  waves  with  microrotation 
(longitudinal  or  transverse  waves).  As  far  as  I  am  aware,  these  waves  have  not 
yet  been  demonstrated  in  the  ultrasonic  region.  Special  transducers  will  have  to 
be  designed  to  show  this  effect. 

4 . 8  Shock  waves 

These  are  the  ones  most  used  for  composite  materials.  Various  methods  have 
been  adopted,  such  as  those  of  Whittier  and  Peck*^,  Sve  and  Okubo*^,  Tauchert  and 
Moon^.  The  same  techniques  are  used  as  for  conventional  solids.  Amongst  others, 
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these  include  methods  using  a  mechanical  projectile,  explosive  charges,  shock 
waves  in  a  gas  and  waves  in  a  shock  tube.  The  shock  wave  is  propagated  in  the 
composite  material  and  detection  of  the  wave  is  done  on  the  opposite  face  by  a 
capacitative  displacement  pick-off,  Fig  18.  Wire  resistance  strain  gauges  can  be 
embedded  in  the  composite  material  to  record  the  waves  at  various  points. 

Large  amplitude  transitory  waves  can  be  generated  by  the  'flying  plate 
projectile'  technique.  This  plate  is  made  from  a  homogeneous  material,  mechanical 
and  polished  to  form  the  head  of  the  projectile,  which  is  launched  by  a  gas  gun. 
The  composite,  in  the  form  of  a  thick  block,  is  attached  to  two  adaptors 
(’buffers')  which  serve  as  connectors  and  which  have  been  aligned,  previously, 
with  respect  to  the  flying  plate.  Displacement  measurements  can  be  by  the  short- 
circuit  contact  method. 

Piezo-electric  crystals  can  be  used  or  a  laser  interferometer  (Michelson, 
with  Doppler  effect)  which  measures  the  displacement  through  a  window*^ . 

Amongst  the  theories  given  earlier,  it  is  the  theory  of  media  with  inter¬ 
action  which  is  seen  to  be  the  most  successful  in  the  study  of  shock  waves 
(equation  (18)).  The  slower  progress  with  the  other  theories  is  explained  by  the 
fact  that  these  latter  contain  parameters  from  undetermined  behavioural  laws. 

The  evaluation  of  these  constants  (or  functions)  by  analytical  or  experimental 
methods  is  very  critical,  although  certain  attempts  (such  as  those  mentioned  in 
section  4.1)  are  promising.  Fig  19  shows  the  good  agreement  between  theory  and 
experiment  in  a  shock  test  on  carbon-phenolic,  resin  composite  materials. 

4 . 9  Mechanical  tests  on  bones 

Bones  are  really  biological  composite  materials.  Experimental  and  theor¬ 
etical  studies  devoted  to  composite  materials  can  be  transposed  to  the 
biomechanics  of  solids. 

The  anisotropy  of  bones  is  very  pronounced  and,  moreover,  is  often  compli- 
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cated  by  a  range  of  elastic  constants.  Guynemer  has  made  measurements  to  give 
the  mechanical  characteristics  of  bones. 

Within  this  field,  subsequent  special  studies  should  take  into  account  the 
specific  characteristics  of  bones.  Rayleigh  waves  in  cylinders  could  be  used  to 
study  cortical  layers. 

5  CONCLUSIONS 

Taking  into  account  the  fractional  aspects  of  the  composite  materials 
considered  here,  it  is  difficult  to  draw  any  firm  conclusions.  At  best,  it  can 
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be  said  that  in  theory  numerous  formulae  can  be  found  for  the  study  of  waves  and 
vibrations.  The  most  simpLe  versions  of  the  theory  of  elasticity  in  microstruc¬ 
tures,  known  as  the  effective  stiffness  theory,  are  capable  of  being  used  with 
success  in  multilayer  composite  materials.  Variational  methods  are  very 
promising.  The  same  is  true  of  asymptotic  methods  and  the  theory  of  media  with 
interaction. 

To  my  mind,  the  variety  of  theories  available  could  cover  all  the  experi¬ 
mental  requirements  without  the  necessity  for  further  development.  Any  lack  is 
with  the  experimenters,  because  the  absence  of  exact  solutions  in  numerous 
problems  requires  experimentation  to  decide  between  the  approximate  theoretical 
solutions  proposed.  It  is  likely  that  more  and  more  they  will  be  forced  to  use 
theories  with  microstructures.  The  use  of  multiple  scales,  one  over  all  and  the 
other  for  the  elementary  lattice,  is  an  approach  which  compels  recognition,  more 
and  more,  to  show  up  certain  fracture  or  delamination  phenomena  which  are  of  a 
microstructural  order.  To  this  end,  attempts  to  remain  within  the  framework  of 
the  classical  theory  of  continuous  media  are  already  reaching  their  limitations. 
A  more  careful  study  of  the  field  of  constraints  by  means  of  a  microstructure 
theory,  of  higher  order  than  the  existing  ones,  should  assist  experimenters  to 
give  a  better  explanation  of  the  phenomena  which  arise  at  the  matrix-fibre 
interface  level. 

Various  physical  methods  (in  particular  optical  methods  and  those  using 

very  high  frequency  ultrasonics)  have  not  been  used  for  composite  materials. 

These  new  methods,  as  yet  little  used  by  engineers,  could  be  capable  of  explain- 
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ing  certain  acoustic  phenomena  already  predicted,  theoretically,  by  Eringen 

Hence  composite  materials  offer  a  wide  field  of  choice  for  interdisciplin¬ 
ary  work  between  physicists  and  engineers. 

Finally,  at  the  calculation  level,  the  question  is  to  know  whether  or  not 
the  various  existing  theories  can  produce  numerical  solutions  which  are  more 
etiective  than  those  obtained  by  direct  examination  of  material  in  finite 
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Table  1 

SUMMARY  OF  THE  ASSESSMENT  CALCULATION  METHODS,  TAKEN  FROM  CHEVALIER 


Nature 
of  the 
method 

Author 

Nature  of 
the  composite 
material 

Advantages 

Disadvantages 

Variational 

Hill 

Isotropic . 
Slightly 
different 
charac  t eristics 

Limits  indepen¬ 
dent  of  material 
geometry.  Docu¬ 
mentary  results. 
Any  interfacial 
bonds.  Pedago¬ 
gical  interest 

Limits  too  wide 
if  the  phases 
are  different 

(limits  for 
modules) 

Hashin- 

Rosen 

Transverse 

isotropic 

Use  of  material 
geometry.  Docu¬ 
mentary  results 

Relatively 

complicated 

calculations 

Hashin- 

Shtrikman 

Isotropic  (and 
even  multi¬ 
phase) 

Limits  closer 
to  those  given 
by  Hill 

Demonstration 
difficult, 
needing  calculus 
of  variances 

Pickett, 

Leissa 

and 

Clausen 

Tetragonal 
orthotropic  or 
transverse 
isotropic 

Applicable  to  a 
large  number  of 
materials. 
Rigorous  method 

Needs  a  large 
computer.  Numer¬ 
ical  results. 

Long  calcula¬ 
tions,  errors 
possible 

static 

Chevalier- 

Vinh 

Transverse 

isotropic 

Very  simple 
documentary 
resul ts . 

Possible  exten¬ 
sion  to  the 
viscoelastic 
region 

Approximation  of 
the  limiting 
conditions.  Only 
applicable  to 
transverse 
isotropic 
materials 

Elasto- 

dynamic 

Behrens 

Anything,  but 
but  must  be 
periodic 

Very  simple 
documentary 
results . 
Applicable  to 
a  large  number 
of  materials 

Transposition  to 
the  viscoelastic 
region  difficult. 
Inapplicable  to 
short  wave  case 
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Table  2 


DIFFERENT  VALUES  OF  THE  ELASTIC  MODULI  OF  A 


SILICON-PHENOLIC  RESIN  COMPOSITE  MATERIALS 


WITH  63%  OF  FIBRES  (UNITS  daN/mm  ) 

Cijk!l 

Chevalier- 
Vinh*  ^ 

Pickett 

solution 

Behrens^ 

solution 

Ultra-sonic 

measurements 

solution 

Vinh 

C1  1  1  1 

1616 

1684 

1651 

1740  ±  30 

C1  122 

718 

667 

683 

756  ±  50 

C 1  1  32 

569 

571 

569 

680  ±  80 

C3333 

5140 

5141 

5142 

4200  ±  100 

C2323 

515 

519 

515 

570  ±  12 

^1212 

449 

509 

483 

492  ±  10 

Fibre: 

Matrix: 


7500  daN/mm 
371  daN/mm 


! 


Table  3 

SOME  POSSIBLE  METHODS  FOR  PHYSICAL  MEASUREMENTS 


Physical  methods 


Remarks 


X  rays 


Neutron  diffraction 


Measurements  in  K. . 

ij 


Measurement  of  specific  heat 


Infrared  spectroscopy 


Transparent  bodies 

Measurements  of  K. . 

ij 


Raman  spectroscopy 


Transparent  bodies 

Measurements  of  K. . 

ij 


Ultra-  and  hyper-sonics 


Measurements  of  C...  „  and  K. .  coupling  terms. 

1  jkf.  lj 

Coupling  between  translation  and  molecular  modes 


Brillouin  diffusion 


Crystals.  Faraday  effect,  rotational  ability  (?) 
Coupling  between  the  subdued  optical  mode  and 
acoustic  phenomena 


For  anisotropic,  micropolar  media  the  elastic  constants  are  defined  by 
the  expression  for  deformation  energy: 


W  = 


|C..,  e..e  +  jK..(r.  -  iJj.)  +  jf..  e..(r 
1  j  k  ?.  ljkS.  1  j  i  Y  l  ljkij  k 


V 


where  e..  =  infinitesimal  deformation 

lj 


r . 

=  |e.  (u  -  u  ) 

=  macrorotation 

i 

lmn  n,m  m,n 

h 

=  microrotation 

e . 

=  permutation  tensor 

iran 


permutation  tensor 
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Figs  1&2 


Key: 

Matdriaux  Composites 
Gdologie  (etc) 

Physique  des  Mdtaux  (etc) 
Plasturgie  Chimie,  etc 
Optique  (etc) 

Acoustique  (etc) 
Biomdcanique  (etc) 

Physique  Math£matique  (etc) 
Mdcanique  (etc) 

Mdcanique  Quantique  (etc) 


=  composite  materials 

=  geology  (propagation  through  stratified  media) 

=  metal  physics  (new  materials) 

=  plastics.  Macromolecular  chemistry 
=  optics  (MoirS-holography,  etc) 

=  acoustics  (vibrations,  ultrasonics,  hypersonics) 
=  biomechanics  (bones) 

=  mathematical  physics  (disturbances,  etc) 

=  mechanics  (micropolar  continuous  media) 

=  quantum  mechanics  (network  dynamics) 


Fig  1  Outline  of  interdisciplinary  studies  which  affect  composite  materials 


a)  Repartition  hexagondle 

Key: 

Ma  trice  =  matrix 
Fibre  =  fibre 

Vide  =  empty  (space) 


b)  Repartition  aleatoire 


Fig  2  Cross-section  of  a  composite  material  in  the  plane  normal  to 
the  fibre  axes 

(a)  hexagonal  distribution 

(b)  random  distribution 


4 


Figs  3&4 


Flolules  C1J  en  daN/mm2  =  C1:|kl  modules  in  daN/mm2 
Concentration  de  fibres  S2>  %  =  concentration  of  s^  fibres,  % 

-  =  Cheval  ier-Vinh15  solution 

o  =  Pickett  numerical  solution 

-  =  Behrens16  solution 

Fig  3  Elastic  constants  C1^  of  a  composite  material  with 
hexagonal  symmetry 

Fibre:  E£  =  7500  daN/mm^,  V2  =  0.2 
Matrix:  E  =  371  daN/mm2,  v-|  =  0.34 


Ke^: 

Terre  dpoxy  =  glass  epoxy 

o  =  experimental  points 

-  =  upper  limit 

-  =  lower  limit 

Fig  4  Determination  of  the  limiting  surface  of  a  multi-layer  glass-epoxy 
composite  material  by  limit  analysis 

The  index  f  refers  to  the  fibre  from  Le  Nizerhy23 

The  errors  between  theory  and  experiment  are  due  to  manufacturing 
imnerfections  and  to  the  influence  of  interfacial  adherance 
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Figs  5&6 


Key: 

Direction  de  propagation 

Plan  de  1'onde  longitudinale  (etc) 

Plan  de  1'onde  a  microrotation  (etc) 

Onde  de  ddplacement  (^3) 

Onde  a  microrotation  (v^) 


=  direction  of  propagation 
=  plan  of  the  longitudinal  wave  (speed  v-, ) 
=  plan  of  the  wave  with  microrotation 
(speed  v£) 

=  displacement  wave  (v^) 

=  wave  with  microrotation  (v^) 


Fig  5  (a)  Waves  with  longitudinal  displacement  and  waves  with  microrotation 
(b)  Coupled  transverse  waves  (speeds  and  v^)  from  Eringen38 


Fig  6  (a)  Curve  of  the  square  of  the  speed  (wave  with  microrotation) 
as  a  function  of  pulse  rate; 

(b)  Curve  of  the  squares  of  the  speeds  and  vjj  for  transverse 

waves  and  waves  with  microrotation  as  a  function  of  pulse  rate  u  . 

38 

The  co  and  o'  pulses  can  be  seen.  From  Eringen  . 
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Figs  7&8 


Fig  8 


Fig  7  Fibre- reinforced  composite  material. 
Overall  and  local  coordinates 


Nombre  d'onde  =  wave  number 

-  =  theory  of  elasticity 

-  =  effective  stiffness  theory 

Comparison  between  the  theory  of  elasticity  and  the  theory 
of  effective  stiffness.  Longitudinal  .wave  in  the  direction 
of  the  layers.  From  Sun  and  Achenbach39  y  = 
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Nombre  d'onde  K 


Fig  9 


Key: 

Vitesse  de  phase  =  phase  velocity 
Nombre  d'onde  K  =  wave  number 


—  =  elastic  solution 

—  =  1st  order  effective  stiffness  theory 

—  =  2nd  order  stiffness  theory 


(a)  System  of  coordinates  at  the  boundaries  for  a  multi-layer 
composite  material 

(b)  Comparison  between  the  three  theories  of  elasticity  and 
1st  and  2nd  order  of  effective  stiffness. 

From  Drumheller^  ;  y  =  u-/u„ 


Figs  iO&Tl 


Key: 

Vitesse  de  phase  =  phase  velocity  (cm/ps) 
Frequence  =  frequency  (Hz) 


-  =  theory  of  mixtures 

-  =  exact  solution 

o  =  experiment 


Fig  10 


Phase  velocity  as  a  function  of  frequency, 
composite  material.  From  Drumheller^' 


Carbon-phenolic 


Key: 

"("see  Fig  10) 


Fig  11  Phase  velocity  as  a  function  of  frequency, 
composite  material.  From  Drumheller*' 


Boron-phenolic 
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Fig  12 


Key: 

VoTtmfetre 
Oscilloscope 
Ponts  de  mesure 
Capteur 
£prouvette 
Pot  vibrant 
Ampl i ficateur 
GSndrateur 
Barre  d'inertie 


=  voltmeter 
=  oscilloscope 
=  measuring  bridges 
=  pick-off 
=  sample 
=  vibrating  pot 
=  ampl ifier 
=  generator 
=  inertia  bar 


Fig  12  (a)  Sketch  of  the  principle  of  an  elasticimeter  with  alternate 

bending,  for  the  determination  of  the  complex  Young's 
modulus  of  composite  materials. 

(b)  The  sample  is  excited  into  alternate  bending  by  a  vibrating 
pot.  The  torsion  bars  at  the  ends  facilitate  the  display  of 
the  vibrations.  Determination  of  Coulomb  moduli.  The 
apparatus  is  similar  to  that  shown  in  Fig  12(a) 


Fig  13 


Key: 

Goniomfetre  =  goniometer 

Eau  =  water 


Rdcepteur  =  receiver 

Oscilloscope  =  oscilloscope 

fmetteur  =  transmitter 

fprouvette  =  sample 

Gdndrateur  d' impulsions  =  impulse  generator 

Ampl ificateur  =  amplifier 

Synchronisation  =  synchronisation 


Fig  13  Principles  of  measuring  elastic  constants  by  means  of  progressive 
ultrasonic  waves 

(a)  Immersion  apparatus.  The  sample  is  fixed  to  a  goniometer. 
Either  longitudinal  waves  or  decoupled  transverse  waves 
are  obtained. 

(b)  Direct  contact  apparatus.  The  decoupled  L  or  T  waves  are 
obtained  from  transducers.  Measurement  of  transit  times 
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Figs  16&1 


Key: 

N ombre  d'onde  =  wave  number 


Transverse  wave: 

-  theory 

•  experiment 


Fig  16  Comparison  between  the  theoretical  and  experimental 
dispersions.  From  Robinson6'.  The  pass-band  regions 
can  be  seen 


Sort  ie 


Frequence  ,  M  Hz 


Key: 

Entrde  =  input 
Frequence  =  frequency,  MHz 
Sortie  =  output 

Fig  17  Fourier  transform  spectral  analysis  of  the  input  and  output 
pulses,  from  Robinson6' 


■ 


I 


Figs  18&19 


Key: 

Adapts teur  d'entrde  =  input  adapter 
Adaptateur  de  sortie  =  output  adapter 
Fenetre  =  window 

Projectile  =  projectile 

Composite  =  composite  material 


Fig  18  Displacement  measurement  (in  a  shock  test)  by 
interferometry 


_ . _ Theorie  des  melanges 


Experience 


Key: 

Vitesse  moyenne  =  mean  velocity  (V/Vq) 
Temps  =  time  (us) 

-  theory  of  mixtures 

-  experiment 


Fig  19  Shock  tests  on  a  carbon-phenolic  composite  material . 
Comparison  with  the  Hegemier  and  NayfehAA»45  theory 
of  mixtures 
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